Fermi edge singularity and finite frequency spectral features in a semi-infinite ID wire 
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We theoretically study a charge qubit interacting with electrons in a semi-infinite ID wire. The 
system displays the physics of the Fermi edge singularity. Our results generalize known results 
for the Fermi-edge system to the regime where excitations induced by the qubit can resolve the 
spatial structure of the scattering region. We find resonant features in the qubit tunneling rate as 
a function of the qubit level splitting. They occur at integer multiples of hvp/l- Here vf is the 
Fermi velocity of the electrons in the wire, and I is the distance from the tip of the wire to the point 
where it interacts with the qubit. These features are due to a single coherent charge fiuctuation in 
the electron gas, with a half-wavelength that fits into I an integer number of times. As the coupling 
between the qubit and the wire is increased, the resonances are washed out. This is a clear signature 
of the increasingly violent Fermi-sea shake-up that accompanies strong coupling. 

PACS numbers: 73.40.Gk, 72.10.Fk 



I. INTRODUCTION 

Systems in which a localized impurity, with an internal 
quantum mechanical degree of freedom, interacts with an 
electron gas, play an important role in many-body the- 
ory. On the one hand, they allow theorists to investigate 
interaction effects and many-body correlations beyond 
the perturbative regime. On the other hand they explain 
observed phenomena such as the resistance minimum (as 
a function of temperature) in dilute magnetic alloys, i.e. 
the Kondo effect fi Another impurity phenomenon that 
has been studied extensively is the so called Fermi edge 
singularityi^i^ In its original incarnation, the effect refers 
to power law singularities in the soft-x-ray absorption, 
emission, and photoemmision spectra of metallic sam- 
ples. As with the Kondo effect^ the phenomenon has 
received renewed attention due to technological break- 
throughs in nano-physics and quantum transport. The 
same physics that is behind the Fermi edge singular- 
ity describes for instance tunneling into and out of a 
small quantum dot coupled to an electron reservoir i^ Re- 
cent studies have focused on non-equilibrium^^— non- 
stationary^ and band structure^^ phenomena. 

In this paper, we confine our attention to the equilib- 
rium situation but consider a setup in which the spectral 
function has interesting features at energies away from 
the singularity. The setup can be realized with current 
technology. The system we study consists of an electron 
gas interacting with a two level system (charge qubit). A 
concrete realization of the qubit could be an electron that 
may occupy the lowest two states of a double quantum 
^Q^ 15,16 pi-Q]-Q ^}j(; point of view of the electron gas, the 
qubit acts as a dynamic localized impurity, while from the 
point of view of the qubit, the electron gas acts as a dis- 
sipative environment. The qubit state-space is spanned 
by the vectors |-|-) and |— ). The total Hamiltonian for 
the system is H = Hq + Ht where 



The energy e represents a gate voltage that controls the 
qubit level splitting and 7 is a small tunneling amplitude 
between the two qubit states. Both these parameters are 
typically under experimental control. 

The Hamiltonians H± = T -f V± describe the electron 
gas. The kinetic term T is the same for both Hamil- 
tonians. The finite range potentials V± represent the 
electrostatic potential produced by the qubit. This po- 
tential depends on the internal state of the qubit, so that 
V+ ^V-. We will specify the system in more detail in 

secnn 

The quantity of interest in this article is the qubit 
tunneling rate W , which is defined as follows: Consider 
the situation where initially the tunneling amplitude 7 is 
zero, and the qubit is prepared in the state |+). The elec- 
tron gas is allowed to equilibrate. We assume zero tem- 
perature so that it equilibrates to the Fermi-sea ground 
state 1-^+) of the Hamiltonian H+. At time t = 0, the 
tunneling amplitude 7 is then switched on, and the state 
|-|-) acquires a finite lifetime. Provided that this lifetime 
is long enough that we can still speak of well-defined 
qubit levels, the probability ri+(i) to find the qubit in 
the state |-|-) decays exponentially^^ 
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(1.2) 



where W is the qubit tunneling rate. 

In order to formulate a quantitative criterion for when 
exponential decay occurs, it is useful to define an energy 
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Here E'q are the ground state energies associated with 
the Fermi sea ground states \F±) of H±. For the system 
described by Hq , the minimum energy difference between 
a configuration (of qubit plus Fermi gas) with the qubit in 
the state |-|-) and one with the qubit in the state |— ) is cj. 
Due to the finite lifetime W~^, there is an uncertainty of 
order W in the qubit level splitting. For the energy levels 
associated with |±) to remain well-defined, the minimum 
energy difference uj between a |+) and a |— ) configuration 



must be much larger than this uncertainty, i.e. 

In this regime n^{t) shows exponential decay. 
By applying Fermi's Golden, one obtains 



(1.4) 



W=\jf dte''*{F+\e'"^'e-'"-*\F+). (1.5) 

J —oo 

(See Appendix \X\ for details). We will use this result to 
obtain an explicit expression for W in terms of the energy 
£ and the potentials V± . In the original incarnation of the 
problem, this quantity corresponds to the photoemmision 
spectrum, i.e. the intensity of electrons ejected from the 
metal, at a fixed x-ray frequency e. (See for instance Sec. 
IVof Ref. IS.) 

In the language of the Fermi edge singularity, the quan- 
tity {F+\e'"+*e-'^-*\F+), that appears in Eq.O is 
known as the closed loop factor. (See for instance Sec. 
Ill D of Ref. ll^-) It is known that, for t much larger than 
the time an individual electron spends in the scattering 
region 
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iH+t-iH- 
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where AE — Eq — Eq is the difference between the 
ground state energies of iJ_ and H+ and A is an ultra- 
violet energy scale. (The branch with arg(iAi) — ±7r/2 
is implied.) The power law exponent a is determined by 
the single particle scattering matrices S± associated with 
the fermion Hamiltonians. Explicitly^i^ 



a = Tr 



'inS'+S'i 
2^ 



(1.7) 



For sufficiently small w, the asymptotic form of Eq. ll.6l 
gives rise to a tunneling rate 



W 



r(a) 



(I)°U^M, 



(1.8) 



with 9(uj) the unit step function, i.e. 9{uj) — 1 for uj > 
and 6'(aj) = for a; < 0. This power law remains 
valid while to <C inm{vF/l, Ep, D — Ep}- Here vf is the 
Fermi velocity, Ep is the Fermi energy measured from 
the bottom of the conduction band, and D is the band 
width. The length scale I is the size of the scattering 
region. This is not necessarily the same length scale as 
the range of the qubit interaction potentials V±, which 
we denote a. Consider for instance a qubit placed next 
to a semi-infinite wire. (See Fig.[T]) Here the size of 
the scattering region is the distance from the tip of the 
wire to the point closest to the qubit, which can be much 
larger that the range of the potential produced by the 
charge of the qubit. In general I > a. 

The condition u> <C vp/l comes about as follows. An 
energy uj corresponds to density fluctuations in the elec- 
tron gas with wavelengths at least vpjuj. The result 



of Eq. ll.8l breaks down as soon as this wavelength is 
short enough for these excitations to resolve the spa- 
tial structure of the scattering region. The restriction 
ijj <^ min{i?i?, D — Ep\ is due to the fact that Eq. ll.8l 
becomes invalid when particle or hole excitations are cre- 
ated in the electron gas close to the band edges. 

An analytical expression for A was obtained by Tan- 
abe and Othaka^^ in the regime Ep ^ vpjl^ where the 
wavelength of electrons near the Fermi level are too long 
to resolve the spatial structure of the scattering region, 
so that the potentials V± may be approximated as b- 
functions. (This is referred to as the limit of contact 
potentials.) A was found to be of the order of Ep. In the 
same limit, approximate results for the finite w behavior 
of W has been obtained. As a function of w, these results 
contain features on the scale of Ep that are associated 
with the band structure of the model. For more detail 
the reader is referred to the review [I^l . 

The regime of Ep <^ vp/l applies in a semi-conductor, 
where it is not uncommon for the Fermi wave-length to 
be large compared other length scales in the problem. 
However, the opposite regime, where Ep 3> vp/l, also 
has physical relevance: In a metallic sample, the Fermi 
wavelength is comparable to the lattice constant, while 
all length-scales associated with the potential are much 
larger. We are not aware of any work in which this regime 
is investigated. 

In this article we study the regime where Ep 3> 
vp/a > Vp/l. There are two significant differences be- 
tween this regime and the previously studied regime. 
Firstly, the ultra-violet energy A is no longer of order 
Ep, but rather is determined by the potential vp. Sec- 
ondly, the rate W as function of ui starts deviating from 
the power law of Eq. ll.8|, at energies ^ vp/l rather than 
at energies ~ Ep. The source of the deviations is no 
longer related to band structure, but to excitations re- 
solving the spatial structure of the scattering region. We 
confine our attention to the case of an electron gas in a 
single chiral channel at zero-temperature. We pay par- 
ticular attention to the example mentioned above of a 
qubit interacting with a semi-infinite wire, where I 3> a. 
We were able to obtain exact analytical expressions for A 
and for the closed loop factor at arbitrary times. We find 
that A < vp/a, and that A depends only on the shape, 
not the magnitude, of V±, i.e. scaling V± -^ cV± leaves A 
unchanged. We were able to compute the tunneling rate 
W{uj) away from the threshold cu —>■ 0~^. We find that 
W{uj) has resonant features at an energy scale vp/l, that 
reveal the nature of many-body correlations induced by 
the qubit. Our main results are contained in Eqs. 14.201 
SSSlEl and[5TUl 

Our analysis is based on the approximation of taking 
the limit Ep — >■ (X) and linearizing the electrons' disper- 
sion relation around the Fermi level, while still taking 
into account the full spatial dependence of the poten- 
tials V±. As pointed out by Gutman et al. (footnote 
36 of Ref.flOl) special care must be taken when lineariz- 
ing the dispersion relation in order to account for the 



anomalous contribution to (i^+| e 



iH+t^-iH-t 



\F+) that 



is related to the Schwinger anomalyj^ In the derivation 
that we present, this anomalous contribution appears 
quite naturally.^^ Our results are obtained by means of 
bosonization,^^'^'^ the application of which to the Fermi 
edge singularity was pioneered by Schotte and Schotte.— 
Bosonization maps the Hamiltonian of Eg. 11.11 onto an 
equivalent spin-boson model^ where a spin is coupled 
to a bosonic bath. For the example of a semi-infinite 
wire interacting with a qubit at a point on the wire that 
is a distance I ^ a from the tip of the wire, the bosonic 
bath spectrum has non-trivial structure. This in turn 
is what leads to the non-trivial finite w behavior of the 
tunneling rate W{uj). 

The rest of this article is structured as follows. In 
Sec. In] we specialize to a Fermi gas consisting of a single 
chiral channel, and introduce a model to describe a semi- 
infinite wire interacting with a qubit at a point a distance 
/ from the tip of the wire. In Sec. lIIII we collect the re- 
sults from the theory bosonization that are required for 
our analysis. We also discuss Anderson's orthogonality 
catastrophe from the point of view provided by bosoniza- 
tion. In Sec. II VI we give a general and exact formula for 
the closed loop factor. This allows us to calculate the 
ultraviolet energy scale A exactly. In Sec.|V] we apply 
the general results of Sec. lIVI to the specific system intro- 
duced in Sec. nil for which the tunneling rate W{ijj) has 
non-trivial features at finite lo. 



II. A SINGLE CHIRAL CHANNEL 

Here we give a mathematical definition of the type of 
electron gas we study. Associated with the electrons in a 
chiral channel of length L with periodic boundary condi- 
tions are creation and annihilation operators ip^{x) and 
'ip{x) that respectively create or annihilate a fermion in 
the state \x) localized at position x. They obey the usual 
anti-commutation relations 



{xl>{x),i,\x')} 



E 

n— — oo 



5{x 



nL), 



(2.1) 



and are periodic with period L. At the point in our 
derivation where it becomes convenient to do so, we send 
the system size L to infinity. 

The non-interacting many-fermion Hamiltonians H± 
have the same linear dispersion but different external po- 
tentials. Without loss of generality (see Appendix [B] for 
details), we can set the external potential in H^ to zero, 
so that 

rL/2 
H+= dxij\x){~id^-^i)4>{x), (2.2) 

J-L/2 

while H^ = H+ + V with 

L/2 

V=l dxv{x)p{x), p(x) = V^(a;)V'(a;) (2.3) 

-L/2 



(Here we work in units where the Fermi velocity vp = 
1.) Associated with i/_ is the one-dimensional scattering 
matrix e"™" where 



rL/2 
Vq = dxv(x). 

J-L/2 

The ground state of -ff+ is the Fermi-sea 

i^+>=n4io), 

A;</.t 



(2.4) 



(2.5) 



where |0) is the state with no particles, the operator 

rL/2 
l-L/ 



4 = —/ dxe^^^i^Hx). 
yL J-L/2 



(2.6) 



creates a fermion in a momentum eigenstate and k is 
quantized in integer multiples of 27r/L. 




FIG. 1: The system described by the potential v{x) of Eg . 12.71 

(a) shows the system with the qubit in the initial state |-|-) and 

(b) shows the system after the transition. The shaded rectan- 
gle represents the semi-infinite wire. The U-shaped contour 
indicates the single chiral channel in which electrons propa- 
gate. The diagram is only schematic. In an actual realization, 
the left and right propagating electrons need not be spatially 
separated, i.e. the U-shape of the contour may be squashed 
into a line. The distance I between the tip of the wire and 
the point on the wire closest to the qubit is indicated. The 
qubit is represented as a double quantum dot with a single 
electron in it. The dashed circle indicates the range of the 
potential through which the qubit and the wire interact. In 
the exploded view, the qubit level spacing e and the tunneling 
amplitude 7 are indicated. While the range of the potential 
and the distance I are of similar size in the figure, we will 
investigate the regime where I is much larger than the range 
of the potential in the text. 



The general results we obtain will be applied to the 
case where the electron gas resides in a semi-infinite ID 
quantum wire. In the limit of a large Fermi energy, this 
system is mapped onto Eqs. ll.l[|2.2l and l2.3l through the 
standard trick of "unfolding" , so that coordinates x and 
—X refer to the same spatial point, but to "different sides 
of the road", i.e. '4)'^{—x) and il}'^{x) create electrons at 
the same position but moving in opposite directions. ^^ 
Thus the potential v{x) is symmetric about a; = 0. The 
system is depicted in Figure [T] As a simple model for 
the interaction between the qubit and the electron gas, 
we will take v{x) = u{x — I) + u{x + I) with 



u{x) 






(2.7) 



Here I is the distance from the tip of the wire to the 
point on the wire nearest to the qubit and a is the range 
of the potential produced by the qubit. This choice of 
u{x) allows us to obtain an exact analytical expression for 
W(uj). In the regime where / 3> a we expect qualitatively 
similar results for any choice of u{x) that is localized to 
a region of length ^ a. 



III. BOSONIZATION 

In this section we collect together the known operator 
bosonization results that are required for our analysis. 
For a tutorial derivation, we refer the reader to Ref.[23l. 
Our notation closely follows Haldane's.— We use these 
results to write H+ and H- in terms of bosonic operators. 
This casts the Hamiltonian H into the form of a spin- 
boson model with a structured environment. We also 
calculate the overlap {F+\ F—), where, as stated below 
Ea. ll.3[ \F±) are the many body ground states of H±, 
which will be relevant when we analyze the tunneling rate 
W in Sec.lWl 

The free fermion Hamiltonian (|2.2p together with the 
Fermi sea ground state \F+) is the starting point for the 
bosonization procedure. Associated with density fluctu- 
ations in the fermion system are operators 



27r 



1/2 



k 



(3.1) 



and aj, q — 2Trn/L, n — 1, 2, 3, . . . that satisfy bosonic 
commutation relations 

[ag,aq']=Q, [aq,al,]:^Sq.q'. (3.2) 

The bosonic annihilation operators aq annihilate the 
Fermi sea \F+), i.e. 



aq\F^ 



0. 



(3.3) 



A central (and non-trivial) result of bosonization is that, 
in terms of the bosonic operators, and for fixed particle 
number 



q>0 



(3.4) 



The fermion density p(x) can be expressed in terms of 
the bosonic operators as 

p{x) ^ N/L +^d. [^{x) + <^t (a;)] , (3.5) 

Zn 

where N = J_',„dxp{x) counts the total number of 
fermions and 



q>0 ^ ^^ 



1/2 



(3.6) 



The if operators satisfy the commutation relations 

[^{x),ipix')] = 0, (3.7) 



[^ix),ipHx')] 
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e i 



■ {x — x')—ri 



.(3.8) 



Using Eqs. l3.5l and 13.61 to express the potential V in 
terms of the bosonic operators, and using expression 13.41 
for 7J_|_, we find for iJ_ 
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where 



L/2 



L/2 



dxv{x)e *' 



(3.9) 
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The Hamiltonian _ff_ is diagonalized by completing the 
square. For this purpose we define new bosonic operators 



bq = ttq 



2n\ 
Tq) 



1/2 



2^' 



(3.11) 



that also obey the standard bosonic commutation rela- 
tions. In terms of these operators the Hamiltonian iJ_ 
reads 



q>0 



9i 



where 
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Substitution of H± from Eqs. l3.4l and 13.121 into the full 
Hamiltonian H of Eg. II. II reveals that the system is 
described by the same Hamiltonian as the spin-boson 
model. (See for instance Ref.[23.) A quantity that plays 
a central role in the spin-boson model, is the bosonic 
environment's spectral function which in our notation is 
given by 



J{q) = TTq 



uq_ 
277 



(3.14) 



(Here we have implicitly taken the L —>■ oo limit.) In 
the context of dissipative quantum mechanics, spectral 
functions J{q) ^ g" for q smaller than some large cut- 
off, play an important role. The case with s = 1 is 
known as an Ohmic environment. We see that a po- 
tential v(x) that is peaked around x = 0, for instance 
v{x) = A/[7r(x^ -|- A^)], produces an Ohmic environment. 
When the spectral function has a more complicated form, 
one talks of a structured bath. A structured bath is ob- 
tained by engineering the potential v{x). As we shall 
show in SeclYl the potential v{x) of Eq. 1 2. 71 produces an 
environment with an interesting structure. 

The ground state energy of H- is i5g ' -I- AE and the 
ground state solves bq \F—) ~ 0, or using the definition 
of bq in terms of aq , 



«J^-> = 



27r 
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(3.15) 



From this follows that the normalized ground state of iJ_ 
is the coherent state 
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/-i/2''^''(^)['^(^)+'^*(^)]/2'^ |F+) (3 16) 



For future reference we note that the overlap 
(F-|-| F—) is easily calculated from Eg. 13. 161 The de- 
tails of the calculation can be found in Appendix [Cl The 
result is 



(F+l F-> H 1^ 



q/2 



(3.17) 



where A is the energy appearing in Ea. ll.6[ and, consis- 
tent with Eq.O (cf. Eq.[! 



a 



(fo/27r) 



(3.18) 



For an explicit formula for A, see Eg. 14. 251 The fact that 
the overlap tends to zero as L~"/^ is known as the orthog- 
onality catastrophe.— The fact that the same ultraviolet 
energy A appears in the closed loop factor and in the or- 
thogonality catastrophe has previously been established 
(for a contact type potential) by Feldkamp and Davis^^, 
and by Tanaka and Othabei^. 



IV. CLOSED LOOP FACTOR 



(The convenience of including the factor exp(iA£'i), with 
AE given by Eq. 13.131 will become apparent below.) 

Having mapped the system under consideration onto a 
spin-boson Hamiltonian, we can simply quote the answer 
from the literature, namely 



P(i) = exp 



1 



'rf.4^(e-*-l) 
q^ \ / 



(4.3) 



with Jici) given by Eg. [3. 141 (See for instance Eqs. 3.35 
and 3.36 of Ref.l25l. but note that P(i) in that work 
refers to a different quantity than the one in Eg. 14. II of 
the present work.) However, we prefer to give a self- 
contained derivation of this result. This derivation goes 
slightly further than simply calculating P(t), namely, it 
produces a bosonic expression for the operator Q(i) that 
is normal ordered, i.e. in which all creation operators are 
to the left of all annihilation operators. This expression 
may in future prove useful for studying non-equilibrium 
effects. Readers prepared to take Eg. 14. 31 as given may 
wish to skip to the paragraph below Eq. l4.21l 

The starting point of our derivation is to consider the 
time derivative of Q, i.e. 



dtQ{t) 



^ie'"+*Ve-'"+ 



'Q{t). 



(4.4) 



Since the Hamiltonian iJ+ is also the momentum opera- 
tor, exp(±JiJ+i) is simply translation by a distance ±i, 
so that 



Thus we find 
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dxv{x)p{x-t) = V{t). (4.5) 
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(4.6) 
x...xe-^-^*o'^*'^(*'), (4.7) 



where t„i — mt/n^ m — Q, . . . , n. 

Let us now consider one of the individual factors in the 
ordered exponent of Eg. 14. 71 Using Eq. l3.5l to relate the 
density operator to the bosonic operators (p and (p^, we 
find 



At pT+At (•L/2 

dt'V{t') = dt' dxv{x) 

Jt J-L/2 

x<^ + ^dt'[^{x^t') + ipHx-t')] 



L 2tt 
AtvoN/L + A{t + At) - A{t) 
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In this section our goal is to calculate the closed loop 
factor 



Pit)^e'^^'{F+\Qit)\F+), 
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where 



m^l 



\x^-^[ip{x-t) + ^\x~t)\, (4.9) 
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Since operators A{t) and A{t + At) com- 
mute to a c-number we have Q-''[M'^+^t)-M'^)] = 

g-[A(r+At),A(r)]/2g-jA(r+At)gJA(r)^ Explicitly 
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We therefore find 

(4.13) 
When substituted back into Ea. l4.7[ this leads to the re- 
sult 
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We can rewrite this as 
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Factor Fi is what Gutman et al.— (see their footnote 36) 
calls the "anomalous" contribution to the closed loop fac- 
tor. An expression involving the determinant of an oper- 
ator acting on single particle Hilbert space often appears 
in the literature^"® in connection with the closed loop 
factor. In Appendix IeI we show that this determinant is 
equal to the expectation value of factor F2 ■ 

Below we consider factors -Fi and F2 separately. From 
Ea. l4.11l we have 



gi[A(t),A(0)]^^^pJ ^ 
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We write factor F2 in boson normal ordered form. This is 
done to facilitate the calculation of the expectation value 
with respect to \F^). 
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dx^^ix^t). (4.18) 
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Explicitly evaluating the commutator in Eg. 14. 171 we find 

^_,[A(t)-A(0)] ^g-.[st(t)-St(0)]g-i[B(t)-i3(0)] 
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Combining this with the result in Eg. 14. 161 for factor Fi, 



we obtain 
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In the expression for P{t), the infinite system size limit 
can straight-forwardly be taken to obtain 



P{t) = exp 
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in agreement with Eq.| 

The large time asymptotics of Eg. 14. 211 can be ex- 
tracted as follows. Firstly we write InP(i) as 



InP(i) 



dq 



lim 



Vq 

27r 


2e 


-iqt 


- 1 






1 


DO 

dq 


/ 

V 


Vq 

2tt 


2 


Vq 

27r 


2 

-1- 


2tt 



-iqt _ -^ 



-qv 



(4.22) 



This expression is then split up into three terms, 
In P{t) = lim,y^o+ (^1 + T2 + Ta), where 

_e"*«*-l 



T^ = 



dq 



Vo_ 
2lT 



T2 ^ - I dq 

/o 9 



27r 



vo_ 
2-K 



n = 



dq 



Q-qiy-it) 



11. 

2-K 



Vq_ 
2-K 



(4.23) 



The integral in T\ is straight forward and leads to T\ — 
— Q;ln(l -I- itjy) where a — {vq/2-kY as in Eg. 13. 181 and 
consistent with Eg. II. 71 In Appendix |D] we show that 
T3 = ©(i"^) and hence vanishes in the large t limit. 
Term T2 can be written as 



To 



-a I dq 

/o 9 



Vq 



1 



(4.24) 



Thus, for large |t|, we obtain P{t) ~ (iAt) " where 

2 ^ " 



A = lim 



-qy 



exp / dq 

y Jo Q 



(4.25) 



This implies that A is determined by the shape of v{x) 
but not by its overall magnitude: The transformation 
v{x) — >■ cv{x) does not affect A. 



V. SEMI-INFINITE WIRE 

We apply the results of the previous section to the case 
of a semi- infinite wire with v{x) given by Eg. 12. 71 so that 



Vq = Xcosiql)e l"^'". 



(5.1) 



Substitution into Ea. l4.21l then yields 



Pit) = c{(i 



-1 x -a/4 




where C 



1 + (l/a) 



a/i 



P[t) ~ (lAt)-" where 

A= |2\/^ 



Thus in the hmit of large t, 



(5.3) 



This last result could also have been obtained using 
Eg. 14. 251 We see that P{t), and thus also the tunneling 
rate W, grows as Z"/^ for I ^ a, 

Some insight into the origin of this result may be ob- 
tained by considering the result in Ea. l3.17l for the overlap 
{F+\ F—), that displays the orthogonality catastrophe. 

For A given by Eq.[521 \{F+\ F-)\'^ « l"^'^- Therefore 
increasing I mitigates the orthogonality catastrophe. The 
tunneling rate W can be written as (cf. Eq. lAlj) 



W 



2vr|7pE^ 



E, 



(+) 



Ei^)\{F+\n~)f 



where En 



(5.4) 



and \n—) are the energies and many-body 
eigenstates of 7f_, the Hamiltonian that describes the 
electrons when the qubit is in state |— ). The results of 
Appendix [C] can be extended to show that not only the 
ground state to ground state overlap, but every over- 
lap |(i^+| '^— )| where \n—) contains a finite number of 
particle hole excitations on top of |i^— ), scales like Z"/^. 
Thus the scaling behavior of W^ ~ l"^^ can be understood 
as a consequence of the mitigation of the orthogonality 
catastrophe. The effect relies on the phase coherence of 
electrons in the section of the wire between |x| = / and 
a; = 0. Hence it is destroyed if the electron phase is 
randomized by impurity scattering or if there is inelas- 
tic scattering. Thus the increase in W with increasing I 
should persist until I exceeds either the elastic or inelas- 
tic mean free path in the wire. Of course the result is 
also only valid as long as w ^ l/l so that excitations cre- 
ated by the qubit do not resolve the spatial structure of 
the potential. Thus, the larger I, the smaller the energy 
window < w ^ l/l in which the enhancement of W 
due to mitigation of the orthogonality catastrophe can 
be observed. 

The tunneling rate W is calculated by expanding the 
third factor in Eg. 15. 21 in a Taylor series in (1 -I- it/2a) 
and Fourier transforming each term separately. Using 
the identities 



dt 



e"^\l+it/r)-' = 



27r {ruY 
iXsj w 



(5.5) 




/A// 



Ico/n 

FIG. 2: The function 1F2 (f ; f , ^; -(li^f) that determines 
the finite u behavior of the tunneling rate H^, as a function 
of w for various couplings a, for the potential v{x) of Eo . 12.71 
The dotted curve corresponds to a = 1/8, the dashed curve 
to a = 1/2, and the solid curve to a = 1. Oscillations with 
period n/l result from the resonant creation of a single particle 
hole pair by the qubit. At larger a the oscillations are washed 
out due to Fermi-sea shake-up, i.e. the excitation by the qubit 
of a multitude of particle hole pairs with a broad distribution 
of energies. 



for s, t > 0, and 



F(s + n) %-r 



T{s) 



m=0 



we obtain 



W{oj) = CWo{uj)Y^ 



(f). 



n\{a)2 



[-{21ujY 



where Wo{uj) is the I = result 



Woiuj) 



27r|7|2 (2aw)" _2 



T{a) Lo 
The factor {a)2n can be rewritten 



•ew). 



«.)..^2-(|)_ 



a + l 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



Substituting this into Eg. 15. 71 we identify the series as 
the Taylor expansion of the hypergeometric function i_F'2, 
yielding one of our main results 

W{io) ^CW,{oj) 1^2 {j; f , ^;-(M2) . (5.10) 

As a; ^ 0+, 1F2 (f ; f , 2±i; -{lujf) tends to 1, so that 
W{uj) indeed has the expected power law singularity for 
small Lo (cf. Eg. ll.8p with A given by Eg. 15. 31 When lo 
becomes of the order 1/Z, excitations in the wire are able 
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to resolve the spatial structure of the potential v{x) and 
W{uj) starts deviating from simple power law behavior. 
In the weak coupling limit, i.e. small a, the hypergeomet- 
ric function 1F2 ( j; ^, ^^; —{luj)^) reduces to cos(?aj)^, 
and W therefore shows oscillations with period n/l as 
a function of oj. These oscillations can be understood 
as being due to the resonant creation of a single particle- 
hole excitation with an energy that satisfies the resonance 
condition uj = irn/l. (This condition maximizes the sin- 
gle particle matrix element (/i| 1^ |p), where h and p refer 
to the single particle orbitals of the hole and the excited 
particle respectively.) 

As shown in Figure [2l the oscillations become damped 
as a is increased. The damping is a signature of a phe- 
nomenon known as Fermi-sea shake-up. At strong cou- 
pling (large a) , rather than creating a single particle hole 
pair, a large number of particle hole pairs are created. 
This corresponds to many charge density excitations (cre- 
ated by the bosonic operators aj) with a broad distribu- 
tion of wavelengths l/q. As a result there is no clear 
resonance any more, and the oscillations in W{u!) are 
washed out. 

It is also instructive to investigating the regime 1// <C 
Lu. Here the asymptotic behavior of the hypergeometric 
function is 



^ , a a a + l ,, ,, 



r(a) 
r(«/2) 



(2luj 



-a/2 



(5.11) 



If we substitute this into the expression fEa. 15.10]) for W, 
assuming / ^ a, we obtain 



Wiu;)[ 



^»i/' r(a/2) UJ 



H- (5.12) 



This is the same rate as would be obtained from a closed 
loop factor 



P{t) 



it 
2a 



(uo/47r)2 



(5.13) 



Such a closed loop factor could also be obtained by cou- 
pling the qubit to two independent chiral channels, where 
the potential the qubit produces in either channel equals 
u{x) of Eq. l2.7l i.e. one of the peaks in the full potential 
v{x) = u{x + I) + u{x — I). This means that, in the sin- 
gle channel semi-infinite wire, at energies w ;:^ 1/Z, the 
potential ^(a; — I) experienced by right-moving electrons 
and the potential u{x -\- I) experienced by left-moving 
electrons, contribute incoherently to the rate W , as if 
left-movers and right-movers belong to separate chan- 
nels. Could this indicate that electrons refiected at the 
tip of the wire have lost all memory of their in-bound 
encounter with the qubit by the time that they again 
reach the qubit on their out-bound journey? Since the 
electrons undergo no relaxation between encounters with 
the qubit, the answer is "no". Rather, what Eg. 15. 121 in- 
dicates, is that processes in which an individual electron 
wave-packet with width l/w <C I is scattered twice, once 



while incident on the tip, and once after being reflected 
at the tip, are rare and make a vanishingly small contri- 
bution to the rate W{uj). 

As stated in the introduction, W corresponds to the ex- 
ponential decay rate of the probability to find the qubit 
in state |-|-), provided that 1 ^ W{ijj)/uj. We conclude 
this section by investigating when this inequality holds. 
For a < 2, W/oj diverges when w — >• 0"'', and the crite- 
rion for exponential decay is violated at small w. From 
the small u asymptotics W{(jj)/uj ^ |7p(a;/A)"/w^, we 
conclude that, for a < 2, exponential decay with rate W 
occurs when 



^»|7||7/Ar/('""V 



(5.14) 



When a > 2, on the other hand, W{uj)/uj no longer di- 
verges, but rather reaches a maximum value of order |a7p 
at w ~ 1/a. Thus, for a > 2, exponential decay with rate 
W occurs for all w > 0, provided that 



1> |a7p 



(5.15) 



The same regime for exponential decay as in Egs. 15.141 
and 15. 151 was identified more rigorously by Legget et ali^ 
in the context of the spin boson Hamiltonian with an 
unstructured Ohmic bath (corresponding to / = in our 
system). (See their Sec. VII. B, and in particular their 
Eq. 7.17a. Note that the quantity that we denote a is 
twice the quantity that they denote a.) The fact that our 
qubit is immersed in a structured bath does not affect 
the result because the bath spectral function J{q) still 
has the same large and small q asyptotics as in the case 
of an Ohmic bath. 



VI. SUMMARY AND CONCLUSIONS 



In 



{F+\e'"+^e^'"-'^ 



this paper we studied the quantity 
\F+), known in the language of 
the Fermi edge singularity as closed loop factor, for 
the case where the Hamiltonians H± describe electrons 
in a single chiral channel. We investigated the regime 
Ep 3> vp/a where a is the typical length scale on which 
the potentials V± associated with H± vary, a regime 
not studied before. We investigated a system where 
the Fourier transform of the closed loop factor gives 
the tunneling rate of a two-level system (charge qubit) 
coupled to the electron gas. Under the assumption of 
linear dispersion (cf. Ea. l2.2p we obtained the exact 
expression for the closed loop factor valid for arbitrary 
V± and arbitrary times. We studied its large time 
asymptotics, and obtained an exact formula for the 
ultraviolet energy A that appears in Eg. II. 61 Unlike in 
the previously studied regime where A ~ Ep, here we 
find A < 1/a <^ Ep. Furthermore it turns out that A is 
determined by the shape, but not the overall magnitude 
of the potentials V±, i.e. scaling V± — ?► cV± leaves A 
invariant. 

We applied our general results to the example a semi- 
infinite wire. The qubit interacts with the wire at a point 



that is a distance I from the tip of the wire. In this 
system we found that the tunneUng rate W could be 
enhanced without increasing either 7 or a. A decreases 
hke Z~^/^ so that W grows hke l"^^. Thus the tunnehng 
rate W becomes larger the further the qubit is from the 
tip of the wire. This effect is due to a mitigation of the 
orthogonality catastrophe. It holds as long as I is less 
than the phase-coherence length of electrons in the wire 
and for level splittings < w < 1/L 

Armed with an expression for the closed loop factor 
that is valid also for small times, we obtain an exact ex- 
pression for W for the semi-infinite wire system. The 
finite w features of W probe the spatial profile of the po- 
tential V± at length scales l/w (in units where vp = 1). 
We study how the w dependence of W changes as the cou- 
pling a between the qubit and the electron gas grows. At 
weak coupling (small a), we find that the rate W oscil- 
lates as a function of lo, and that the period is n/l. This 
is due to the resonant excitation of a single particle-hole 
pair in the Fermi sea. The wavelength of the associated 
charge density fluctuation is tt/lo. The resonance condi- 
tion is that an integer number of half wavelengths fit into 
the part of the wire between the tip and the point where 
the qubit interacts with the wire. At strong coupling 
(large a) on the other hand, many particle hole excita- 
tions are created. This is known as Fermi sea shake-up. 
The corresponding density fluctuations have a broad dis- 
tribution of wavelengths and hence there are no clear 
resonances. This results in the damping of the oscilla- 
tions in W{lli) as a is increased. One of our main results 
(Eq. l5.10l) is an exact formula for this damping by means 
of Fermi sea shake-up. The result is illustrated in Fig. [31 
We also analyzed rate W in the ui ^ l/l limit and saw 
that here the left and right moving electrons contribute 
to the rate W as if they belong to independent chan- 
nels. This happens despite the fact that each electron 
incident on the tip encounters the qubit twice, once be- 
fore being reflected at the tip, and once afterwards, and 
no electron relaxation occurs between qubit encounters. 
The result therefore indicates that processes in which an 
individual electron wave-packet with width 1/w <^ I is 
scattered twice, once while incident on the tip and once 
after being reflected at the tip, make a vanishingly small 
contribution to the rate W{uj). 



Appendix A: Obtaining W from Fermi's golden rule 

In this Appendix we apply Fermi's golden rule to ob- 
tain the expression in Eg. II. 51 for the transition rate W. 
The initial state for the transition is \F+) (g) |-|-) with 
energy E^ -t- e. Possible final states are of the form 
\n—) (E) |— ), where \n—) is an eigenstate of H- and has 
energy En'. We have to sum over all eigenstates of H^. 



Thus 



n 

/oo 
-00 



x\n~)(n~\e-'"-'\F+) 



|7r / dte''' {F+\ e'-«+*e-*"-' \F+) . (Al) 



Appendix B: Gauging away V+. 

In the main text we chose the potential V+ as zero, and 
stated that this does not involve any loss of generality. 
Here we prove this claim. Suppose 

rL/2 
H±^ dxi;\x){-^^^-^l + v±{x))i^{x), (Bl) 

J-L/2 

Now define position dependent phases 

^±{x) = — f dx' v±{x'), 

J-L/2 



(B2) 



and total phase shifts A± = X±{L/2). Then define a new 
set of fermion operators i}}{x) related to V'(a^) by 



il){x) 



— p*[-^+(2:)-A+a;/L] 



V'(x). 



(B3) 



The operators V'^(^') ^iid ^'i^) obey the same anti- 
commutation relations as ijj^x') and ^"(2^) ^-nd are also 
periodic with period L. 

In terms of i'{x) and '(/'^(x), the Hamiltonian iJ+ has 
the form 



H, 



while 



L/2 



dx4!\x){-id^- ljL-X+lL)%l}{x), (B4) 



-L/2 



H~ 



L/2 



-L/2 



dx'ip\x){^ — idx — 11 — X+/L 

+ v-{x)~v+{x))^{x), (B5) 



Thus, in terms of the new fermion operators, iJ+ and 
_ff_ are of the same form as in Eas. l2.2l and 12.31 with 
v{x) — > v^{x) — v+{x) and fi ^>- fi + X+/ L. 



Appendix C: Anderson's orthogonality catastrophe 

Anderson^^ states that the overlap {F+\ F—) vanishes 
as a power law L~"/^ as the system size grows. For the 
present system we can calculate this overlap exactly for 
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arbitrary potentials v{x). Our starting point is Eg. 13. 161 
and the operator identity e^"*"^ = e"^e^e~[^'^l/^, pro- 
vided that [A, [A, B]] = [B, [A, B]] = 0. 



(F+l ^- 



27r 



Lq 



exp ■ 



q>0 

-E 



1/2 



<^+i»pi: S S".-S«: 



2tt 



27r 



27r 



(CI) 



In the large L limit, we rewrite this as 



= lini 



1 

— exp 

y 



dqe 



-qy 



X 2/ exp 




-a/2 



(C2) 



where a = |wo/27rp as in Eg. 13.181 The y — > 0+ Umit of 
the two factors marked 1 and 2 can be taken separately. 
Referring back to Eg. 14. 251 we identify the factor marked 
1 as the energy A that appears in Egs. ll.6l The sum in the 
exponent of the factor marked 2 is the Taylor expansion 
of the logarithm function and hence 



2/exp f ^ —e «M = y [1 - cxp(-27rj//L)]" 



L/27r. 



This leads to the result 

(F+l F-) = 



27r 
AL 



a/2 



(C3) 
(C4) 



Appendix D: Asymptotics of P{i) 

In the main text, in the derivation of the asymptotic 
form of P(i), we stated that T3 in Eg. 14. 231 vanishes like 
\/t in the large |i| limit. Here we give a proof. By Fourier 
transforming from Vq to v{x) we obtain 



T. 



dx i dx 



,v{x) v{x') 



27r 27r 

g-ig(t-iy) 



dq ■ 







q 



iq{x- 



- 1 



The integral / can be performed to obtain 

r , /, «A , f, i{t-x + x') 
/ = In 1 + - - In 1 + -i - 



y/\ y 

Expanding in 1/i we find I = [x — x')/t + ©(i^^) 



(Dl) 



(D2) 



Appendix E: Determinantal formula related to 
closed loop contribution 

Here we show that the expectation value of the factor 
F2 in Eg. l4.15l with respect to |-F-|-) equals a determinant 
of an operator acting on single particle Hilbert space. 

The proof relies on the following general result for 
fermionic systems. Let B = {\n) \n = 1, 2, . . .} be a 
set of orthonormal single particle orbitals and let cj^ and 
c„ be the associated fermionic creation and annihilation 
operators. Let F be a subset of B. Without loss of gen- 
erality, we may take F = {\n) |n = 1, 2, ... N}. Let \F) 
be the many-fermion state 



N 



\F)=Y[cl\0) 

m— 1 

Let H be the operator 

00 
11 — ^ ^ n^jijiCu^c 



(El) 



(E2) 



771, n—l 



Then e*^ |F) — n7n=i ^L |0), where the fermionic oper- 
ator c\^ creates a particle in the orbital \fi) = e*'* |n), 
where 



h^ 



m.n— 1 



hmn \m) {n\ , 



(E3) 



is an operator acting on single particle Hilbert space. 
This imphes that (F| e*^ \F) = det e^pp is a Slater deter- 



minant where e^p is sm N x N matrix with entries 



[efFi 



{m\ e'" \n) 



(E4) 



^N 



The operator / = X]77=i 1*^) ('^l projects onto the sub- 
space spanned by the set F, and hence the matrix e^pp 
has the same determinant as the operator 1 — / -I- e^'Vi 
leading to the result 



(F|e^^|F)=det(l-/ + e"7) 



(E5) 



We derived this result for a state \F) containing a finite 
number of particles. We postulate that a similar result 
holds for the state |i^-l-) representing an infinitely deep 
Fermi sea. 

In order to apply the above result to F2 ~ g-'^l^W-Mo)] 
we have to show that A{t) — A{0) is quadratic in fermion 
creation and annihilation operators. This is indeed so, as 
is seen by referring to Eg. 13. 51 and Eg 14.91 to obtain 



A{t)-A{0) = / dt' 
Jo 



L/2 



dx 



L/2 



VO 



v{x — t') — 'ip{x)^^{x). 



L 



Thus we may use Eq. lE51 to write 

(F+|e-'l^(*)-^(°)l|F-f)=det(l 



-j5. 



(E6) 
(E7) 
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where n and S are the single particle operators 
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